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VORTEX SIMULATION OF THE PRESSURE FIELD OF A JET 

•k 

Y. T. Fung and C. H. Liu 
NASA Langley Research Center 
Hampton, Virginia 

ABSTPvACT 

Fluctuations of the pressure field of a jet are simulated numerically 
by a flow model consisting of axisynmetric vortex rings with viscous cores 
submerged in an inviscid uniform stream. Vortex shedding time intervals, 
randomly created to imitate the time-history characteristics of the pres- 
sure signals of a jet, are generated based on a probability distribution 
of the intervals between successive pressure peaks obtained from experiments. 
It is found that, up to five diameters downstream of the jet exit, the charac- 
teristics of the pressure fluctuations and the most probable time inter- 
vals between experimental and numerical results show good qualitative agree- 
ments. The role played by the axi symmetric vortex model in pressure field 
as well as extensions of the model is also discussed. 


* NAS-NRC/MASA Research Associate 


Symbol s 
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D Diameter of the jet nozzle 

E Complete elliptical integral of the second kind 

Bessel function of the first kind of order zero 
Bessel function of the first kind of order ono 
K Complete elliptical integral of the first kind 

P Point of interest 

p Fluctuating pressure 

p^ Unperturbed pressure 

R Radial position of vortex rings 

r Radial coordinate 

r^ The shortest distance from P to the center of vortical cores 

r2 The longest distance from P to the center of vortical cores 

t|^ Shedding time of the k-th vortex ring 

Convection velocity of vortex rings 

U Velocity of the uniform stream 

Induced velocity in radial direction 

Induced velocity in axial direction 

Z Axial position of vortex rings 

Z Downstream distance for destruction of vortex rings 

max 

z Axial coordinate 

a Decaying rate of the circulation of vortex rings 

r Circulation of vortex rings 

6 Effective radius of vortical cores 

V viscosity 
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p Density of the fluid 

T Time scale 

(}> Velocity potential 

Stream function 
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INTRODUCTION 

The concept of modeling turbulent jets in terms of swne well- 

defined patterns has been proposed by many researchers as a way to 

attack the problem of jet noise generation and its suppression. 

Based on the experimental observations that pressure fluctuations 

outside the mixing region come in rather well -organized packages, 

Mollo-Christensen^ first suggested that certain natural organized 

structures might exist within the shear flow turbulence and that 

these structures might be the primary mechanism for the generation 

of jet noise. A number of subsequent investigations have focused 

on this approach with the hope that this regular pattern, if any, 

and its role in jet noise production could be better defined. 

2 

Later, Becker and Massors carried out a study on the vortex evo- 
lution and instability of an ax i symmetrical jet with a Reynolds 
4 3 

number range of 10 . Crow and Champagne performed measurements 

5 

at a Reynolds number around 10 and reported that a large-scale 
orderly structure was observed within the noise-producing region 
of the jet. Based on these observations, they proposed a theore- 
tical flow model consisting of axi symmetric vortex trains, with a 
preferred Strouhal number of 0.3 emerging by calculation. This 

model of an axial array of vortices was also recommended by other 
4 5 6 

researchers ’ * as the basic structure of a flow pattern designed 
to imitate the relatively periodic and deterministic behavior of 
turbulent round jets. 
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It is the purpose of this numerical investigation to verify 
and analyze the connection between the vortex model and the pressure 
field of a jet. There have been numerous studies of vortex rings; 
however, most of these concentrate only on a single vortex. The 

7 

classical vortex ring of small cross-section in a perfect fluid 
was physically unsatisfactory until the viscous effects were taken 
into account. By considering an inner viscous core and applying 

o 

the boundary -layer technique, Tung and Ting removed the mathema- 
tical singularity which is the essential drawback of the inviscid 
model. The same solution was obtained through an independent ap- 
proach based on considerations of the kinetic energy and momentum 

g 

balance of the fluid . The interaction of two inviscid vortex rings 

was studied by Sommerfeld^^; the viscous case was examined analy- 

12 11 
tically by Gunzburger and experimentally by Fohl and Turner . 

By utilizing the viscous solution to formulate the vortex model 

into a train of axisynwnetric vortex rings with viscous cores, 

13 

Liu et al . developed a combined theoretical and experimental pro- 
gram and obtained good qualitative comparisons between the estima- 
tion and the measured fluctuations of the pressure field of a jet. 
The numerical simulation reported here is undertaken as a further 
test of the validity of the proposed vortex rii^g model. The fluc- 
tuation characteristics of the pressure signals are modeled by the 
random generation of a train of vortex rings. Statistical explana- 
tions to random behaviors of the pressure field of a jet will be 
the main emphasis of the present investigation. 
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Let [R(t), Z(t)] be the position of an Isolated vortex ring 
at time t in an axisytnnwtric cylindrical coordinate system (r,z) 
with z coinciding the axis of symmetry of the ring, and r^ and 
be the shortest and longest distance from the point of interest P 
to the center of the vortical core (see figure 1). Within the 
framev^jrk of the classical inviscid theory, the well-known stream 
fun^tion^ for a potential flow induced by an axisymmetric planar 
vortex ring is 


ip = 





1 


1 


( 1 ) 


Here K and E are the complete elliptical integrals of the first and 
second kind, and r(t) the total circulation around the ring. The 
corresponding induced velocities at P in r- and z-direction are, 
respectively. 


w = 

“l r az 


VI 

•^3 r ar 


( 2 ) 


The well-known deficiencies of this inviscid formulation arise from 
the undefined vortical core and the singularity of the velocity at 
the center of the core. 
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Obtained by considering a viscous inner core in a region of 

g 

large velocity gradients , the equation of motion for a vortex ring 
submerged in an irrotaticnal Inviscid flow stream is 

i = ^ n«9 f - 0.558) (3) 

Here the dot represents the derivative with re.spect to t. The 
symbol 6 is the effective radius of the visa js core defined as: 

5 = 2/ut (4) 

and is assumed to be much smaller than the radius of curvature R. 

The quantity v is the viscosity and t is a time scale defined by 
the integral equation 


t 

x(t) = / R(s)ds/R{t) (5) 

‘o 

where t^ is the creation time of the vortex ring. 

The intention of this numerical study is to sumulate the fluc- 
tuations of the pressure field near a jet by means of a flow model 
consisting of a train of coaxial axisynmetric vortex rings with 
viscous cores submerged in an inviscid uniform stream. The stream 
velocity is either taken as one-half of the jet velocity or 
regarded as a parameter in order to fit the experimental data. 

The measured pressure signals to be simulated, shown in 
figure 2, were obtained from the experiment given in reference 13. 
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The experimental setup (figure 3) consists of a jet nozzle pointing 

upwards and a vertical array having twelve microphones located just 

1 3 

outside the boundary of the jet . In the present simulation, the 
region of Interest 1$ restricted to that of the experimental nwa- 
surements where the convective speeds are small in comparison with 
the sound speed. The Implication Is the Laplace equation d(xn1 nates 
over the acoustics wave equation, and Is, therefore, adequate to 
describe the flow field by the Incompressible potential solutions. 
{T.iis assumption of Incanpressiblllty **s reasonable In the present 
ca:t since compressibility does not alter the turbulence structure 

14, 

until the mean velocity greatly exceeds the sonic speed ). It is 

further hypothesized that the overall contribution to the flow 

field by a train of vortex rings can be approximated by superposing 

the contribution of Individual vortex rings, and hence no merging 

mechanisms will be considered. 

Let the subscripts k and 1 denote the quantities associated 

with the k-th and i-th rings, respectively. As an extension of the 

12 

investigation for the Interaction of a pair of vortex rings , the 
governing equations of motion of the k~th ring in the presence of 
the 1-th one are as follows: 


fi,^(t) = W^^.(t, R,^, Z,^) 


for i f k 


( 6 ) 


N r 8R 

Z^(t) = Wj^Ct, R|^, Z|^) + {log ^ - 0.558} 
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where and defined by equations (2), are the Induced ve* 
loci ties by the presence of other vortex rings. The number of 
active vortices N Is a paranwter to be obtained by matching with 
the experimental data, and Is the uniform stream velocity which 
Is taken as equal to half of the jet efflux velocity, Uj. Assume 
that the k-th rip starts rolling up at the shedding time tj^ with 
a linear strength proportional to U U , U being the convection 
velocity of the ring. The circulation reaches a saturated value at 
the shedding time of the (k+l)-th vortex ring, It Is also as- 

sumed that, beginning at a certain distance downstream of the 
jet, the circulation of the ring decays with a linear rate less than 
the rolling-up strength and that the ring vanishes when its circula- 
tion reaches zero. The basis of this hypothesis for the growth and 
decay of vortex rings will be given in the next section. The evolu- 
tion of the circulation is described by the following equation (see 
figure 4). 


^ f,^(t) 


(7) 


t - t. 


Here f^(t) = < - t^ 


Vi - ^k - 


0 t < t. 


‘k 1 ‘ ^ *k+l 


‘k+1 1 ‘ ' ‘k 


< t < tj* 


tp < t < 
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■ the time whv^n the center of the k-th vortex ring reaches 
the distance downstream. Is a parameter to be 
obtained by matching with the experimental data. 

tjj*- - t|^) + t]J, the time of complete destruction of 

the k-th vortex ring, I.e., * 0. 

a » a parameter which governs the decay of the vortex ring, 
ranging over values between 0 to 1. 

Equations (6) will be solved subject to the Initial conditions 


Rk(tk) = 0/2 

w ■ ° 


(8) 


at the nozzle of the jet. 

The slimjlatlon of the real-time pressure distribution measured 

near a je^ Is to be made along the cone of measurenent In the experl - 
13 

ment , which is outside of the diffusive cores of all the vortex 
rings at all times. The classical Bernoulli solution is then appli- 
cable since the region of interest Is practically irrotatlonal and 
Inviscld. The Induced velocity potential^ for a train of coaxial 
axl symmetric vortex rings Is 


4 , = - 



“ -k!z-ZJ 

/ e sgn(z - Z^) d< 


( 9 ) 



n 


where J and J, are the Bessel function of the first kind of order 0 
0 1 

and 1 respectively. The corresponding pressure variation, obtained 
with the use of the boundary condition that the incompressible uniform 
stream remains unperturbed at a distance sufficiently far away from 
the vortex rings, is given by 


. ifi = 

P 


s + 

3t 


[(u. 


N 

I 

1=1 


w Y 


N 

( S 
1=1 


Wii)'] - 


( 10 ) 


where p Is the density of the fluid (assumed to be constant) and p^ls 
the mean pressure outside the perturbed region. The summations ln 
equatior' (9) and (10) are taken over all the active vortex riros 
w1th1» rj region of Interest. 

The theoretical model under consideration is based on equations 
(6) and (10). The Instantaneous positions of the vortex rings and 
the Induced pressure fluctuations can be obtained by numerical Inte- 
gration of these equations, subject to the Initial conditions (8). 


BASIS FOR THE SIMULATION 

The dominant parameters of the simulation based on the preceding 
mathematical formulation are U^, the uniform stream veloc* y; v, vis- 
cosity; r(t), circulation of the vortex rings; and tj^, the shedding 
time of the vortex rings. These parameters are chosen with reference 
to the experimental observations. 
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Viscosity 

Viscosity enters the simulation in terms of the viscous core 

defined by equation (4). There is not much detailed evidence on the 

growth of vortical cores for the ring type structure of jets except 

that reported by Lau and Fisher®, who estimated the vortex core radius 

to be one-tenth or less of the jet diameter at an axial position two 

diameters downstream of the jet exit. The viscosity, equal or less 
2 

than 0.01 m /sec for the present simulation, is therefore chosen such 
that the vortical core possesses the above-estimated growth rate. 

Circulation 

The time-history of circulation of vortex rings described by 
equation (7) is constructed based on the shear strength in the mixing 
layer close to the jet exit and the experimental observations of the 
vortex ring type structure downstream of turbulent jets. For sim- 
plicity let the convective rate of increase of the circulation equal 

to the average shear strength at the jet exit, i.e., U , and let this 

00 

rate of increase remain constant as the vortex ring travels downstream 
with a convection velocity U^. The circulation reaches a saturated 
value at the shedding time of the successive vortex ring, The 

creation of the circulation described by equation (7) is then followed. 

The assumption on the decay of the circulation of vortex rings 
starting at a certain distance downstream of the jet is based on the 
fact that vortex ring type structures and their corresponding pressure 
signals are seldom observed beyond the potential core. According to 


13 


2 3 

the experimental observation ’ , vortex type patterns start dissipat- 
ing no further than four dian^ters downstream from the jet exit. The 

distance of in this simulation is therefore chosen to be three 
max 

diameters or less. The vorticity then dissipates, under the viscous 
influence, with a linear rate governed by the parameter ot to be ob- 
tained by matching with the experimental data. 

Rando m Shedding Time 

3 

The large-scale vortex puffs were generated randomly in time ; 
an observation which is consistent with the random nature of the meas- 
ured pressure field. The present simulation of the random shedding 
time betwe r successive vortex rings are estimated based on the time 

intervals between successive peaks of pressure of fluctuations as 

13 

measured in the experiment . The simulated pressure is then com- 
pared with the measured pressure fluctuations from which this numerical 
input is obtained. Second, a numerical program is used to generate 
random vortex shedding time intervals having probability distribution 
similar to that of the intervals between successive pressure peaks 
measured in reference 13. This requirement for similar probability 
distribution provides a statistical similarity of peak variations and 
characteristics between the time-history measured pressure and its 
corresponding simulation. 
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RESULTS AND CONCLUDING REMARKS 

Numerical results for the time-history of the pressure in the 
vicinity of a jet, with the time intervals between peaks at one dia- 
meter downstream used as input to the numerical model, are plotted in 
figures 5 along with the experimental pressures. No attempts were 
made on matching their amplitudes because of the random nature of the 
pressure fluctuations. Direct comparisons on pressure peaks are im- 
possible even for two pieces of record obtained from the same set of 
experiment. Statistical analysis appears to be the most sensible com- 
parison. 

Figure 6 shows the simulated pressure fluctuations calculated 
by inputing three hundred random shedding time intervals having pro- 
bability distribution similar to that of the intervals between experi- 
mental pressure peaks at the one diameter downstream station. Statis- 
tical comparisons between these numerical results and the experimental 
13 

measurements are made at several stations along the cone of measure- 
ment. Figures 7 display comparisons of the probability distributions 
of the time intervals between calculated and experimental pressure 
peaks for several values of Z/D. Good qualitative agreement on the 
envelops of the probability distribution and the most probable time 
intervals (interval with maximum number of counts) between pressure 
peaks is observed. The experimental samples outnumber the numerical 
ones by the ratios indicated in figures 7. Improvements on the 
numerical part will be expected if more numerical samples are considered. 
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The variation of the most probable time Intervals between pressure 
peak obtained from the distribution curves in figures 7 Is plotted 
versus the downstream stations along the cone of measurement. The 
numerical results show good agreement with their experimental counter- 
iparts. No comparison beyond the potential core was made 
since vortex ring type structures are rarely seen there. As shown 
by the experimental data, the slope of the most probable time Intervals 
beyond the potential core is different from that within the potential 
core. These two slopes intersect at a distance 4D downstream with a 
corresponding Strouhal number equal to 0.3, an Important location and 

3 

number for the emission of sound . This suggests a possible geometry 
changed for the axisynmetric structures. Nonas 1 symmetric and nonplanar 
C(Mnponents are believed to play an imoortant role beyond the potential 
core, and are now being incorporated into the model in order to allow 
simulation further downstream of the jet exit. 
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Fiqure 2 . Time-history of the measured pressure signals 
at several stations. 
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Fiqure 3. Experimental setup. 
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Figure 5a. Pressure comparison between the simulation and 
the time-history measurements at 3D downstream. 



the time-history measurements at 4D downstream. 



Figure Sc. Pressure comparison between the simulation and 
the time-history measurements at 50 downstream. 



Figure 6. Simulated pressure fluctuations by randomly 
shed vortex rings at downstream stations, 
(U^ =110 m/sec). 
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Figure 7b. Statistical comparisons of the probability distribution 
of the time interval between successive pressure peaks. 




Figure 8. Variations of the most probable peak time 
interval along the cone of measurements. 



